Effective Evolution Equations from 
Many Body Quantum Dynamics 



Benjamin Schlein 

Institute for Applied Mathematics, University of Bonn 
Endenicher Allee 60, 53115 Bonn, Germany 
benj amin . schleinOhcm . uni-bonn . de 

December 22, 2010 
Abstract 

In these notes we review some recent results concerning the derivation of effective 
equations from first principle quantum dynamics. In particular, we discuss the deriva- 
tion of the semi-relativistic Hartree equation for the evolution of boson stars, and the 
derivation of the Gross-Pitaevskii equation for the dynamics of Bose-Einstein conden- 
sates. 

1 Introduction 

Systems of interest in physics are typically composed by a huge number of elementary 
particles. Dilute samples of Bose-Einstein condensates contain 10^ — 10^ atoms (and, strictly 
speaking, each atom contains many elementary components). The number of molecules 
contained in chemical samples is typically of the order of Avogadro's number, Na — 6 • 10^^. 
Systems of relevance in astronomy and cosmology (like stars and galaxies) are composed by 
up to lO^'^ elementary components. 

In principle, the dynamics of these systems can be determined by solving fundamental 
evolution equations like the Newton equation or the many-body Schrodinger equation. In 
practice, however, fundamental equations are impossible to solve (neither analytically nor 
numerically) when so many particles are involved (unless the interaction among the particles 
is neglected). Moreover, observers are not interested in determining the precise evolution of 
every particle. Instead, they need a prediction for the macroscopically measurable properties 
of the dynamics (which result by averaging over the many particles in the system) . For this 
reason, it is very important to find effective evolution equations which, on the one hand, can 
be easily solved (numerically), and, on the other hand, accurately predict the macroscopic 
behavior of the system under consideration. One of the main goal of statistical mechanics 
consists therefore in the development of effective theories approximating the solutions of 
fundamental evolution equations in the relevant regimes. In these notes, we are going to 
discuss two examples of systems of interest in physics, for which the derivation of effective 
evolution equations can be made rigorous in a mathematical sense. In Section [2l we will 
illustrate the derivation of a semi-relativistic Hartree equation for the evolution of boson 
stars. In Section [3l we will sketch the derivation of the Gross-Pitaevskii equation for the 



dynamics of initially trapped Bose-Einstein condensates. In both cases, the starting point 
of our analysis is the fundamental many-body Schrodinger equation for bosonic systems. In 
the rest of the introduction, we show in an abstract setting how effective evolution equations 
emerge from many body quantum dynamics in certain regimes. 

We consider quantum mechanical systems of N spinless bosons in three dimensions (the 
spin does not play any role, and therefore will be neglected). We describe these systems 
on the Hilbert space 71^ = -Lg(M^^), consisting of all functions in L^(R^^) which are 
symmetric with respect to arbitrary permutations of the particles. The time evolution 
is then described by the particle Schrodinger equation 

idttpN,t = HNtpN,t (1) 
for the wave function ^l^N,t S T~{-N- We consider Hamilton operators of the form 

N N 

Hn = Y^ -A,^ + a ^ V{xi - Xj) , (2) 

j=l i<j 

where A G M is a coupling constant and where the potential V{x) describes the (two-body) 
interaction (the precise form of V depends on the system under consideration). Here and in 
the following, we choose units so that Planck's constant h = 1 and the mass of the particles 
m = 1/2. We restrict our attention to (approximately) factorized initial data, where all 
particles are essentially described by the same orbital. The kinetic energy is then of the 
order A'^, while the potential energy is of the order AA^^. To obtain a non-trivial effective 
evolution equation in the limit of large A^, we have to assume that XN^ ~ A^, hence that 
K := AA^ is a quantity of order one. This regime is known as the mean- field regime. 

To analyze the mean-field regime, consider, at time t = 0, the factorized initial data 
il^N,t=o = V''^^; an arbitrary one-particle orbital if G L^(R^). Because of the interaction, 
factorization is not preserved by the time evolution. However, if A^ ^> 1 and k = NX 
is of order one, the interaction is very weak and one may still expect factorization to be 
approximately (and in an appropriate sense) preserved: 

N 

'^N,t^Y{Mxj) (3) 

for an appropriate evolved one-particle orbital ^pt G L^(R^). If this is true, it is very easy 
to derive the self-consistent nonlinear Hartree equation 

idtift = -A^t + K-iV *\ipt\'^)ipt (4) 

for the one-particle orbital ipt- This simple argument suggests that the Hartree equation 
gives an effective description of the evolution of initially factorized bosonic systems in 
the mean-field regime characterized by A^ ^ 1 and fixed k := NX. 

To obtain a mathematical precise statement, we need to specify in which sense ([3]) 
holds true. To this end, we define, for A; = 1, . . . , A^, the /c-particle reduced density matrix 
associated with ipN,t by taking the partial trace 

^nI = Trfc+l,...,7V \^N,t){'>pN,t\ 
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where \il^N,t){ipN,t\ denotes the orthogonal projection onto ipN,t- In other words, 7]^^ is 
defined as a non-negative trace-class operator on L^(M^'^) with kernel given by 

7^](xfc;x'^) = j d-^N~ki'N,t{y^k,y^N-k)i'N,t{^k^^N-k) , 

where = (xi, . . . , Xfc), x'^ = (x'^, . . . , x'^), x^r-fc = (2:^+1, • • • , xn)- Note the normalization 

The next theorem, which holds under suitable assumptions on the potential V , tells us 
that ([3]) can be understood as convergence (in the limit of large N) of the reduced densities 
associated to il)N,t- 

Theorem 1.1. Let if G H'^(M.^), \\ip\\ = 1, n £R and let '4)N,t = e'^^^^Lp®^ be the solution 
of the Schrddinger equation ^j, with initial data ipN = o-nd with Hamilton operator 

N N 

Hn = Y, -^x, + ^ E ^(^^ - ^^■) • (5) 

j=l i<j 

(k) 

Then, ifjj^[ denotes the k-particle reduced density associated with ipN,t, we have, for every 
fixed A: G N and t G M, 

as N ^ 00. Here, the convergence is in the trace norm and (ft is the solution of the Hartree 
equation with ipt=o = 

The first proof of Theorem ll.ll has been obtained by Spohn in [20] for bounded potentials. 
In [8], Erdos and Yau proved Theorem 1 1.1 1 for the (attractive or repulsive) Coulomb potential 
V{x) = ibl/|x|. In P3I, a joint work with I. Rodnianski, we considered again the Coulomb 
interaction, but this time we obtained bounds on the rate of the convergence. In [11] . 
Knowles and Pickl extend the theorem to more singular potentials (with control of the rate 
of convergence). 



2 Dynamics of Boson Stars 

In this section we consider systems of gravitating bosons forming so called boson stars. We 
describe boson stars with the Hamilton operator 

N N 

^g.av = 5:^l-A,,-G5]^-— (6) 

acting on the Hilbert space Hn = ^^(M'^^) (we use a relativistic dispersion for the particles, 
but the interaction is classical). As explained in Section [H we are interested in the dynamics 
generated by ([6]) for large N and small G, with NG of order one. Since the physical value 
of the gravitational constant, in our units, is approximatively given by Gphys — 10"^*^ (for 
bosons with mass comparable to a hydrogen atom), this model can be used to describe 
boson stars with N ~ 10^'' particles. For such values of A^, it makes sense to fix A := NG 
and to study the dynamics generated by 

N . N 

^--E^^-^Er^ (7) 

j = l Kj 
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in the limit N ^ oo. 

The Hamiltonian ([7]) is critical in the following sense. For every N € N, there exists a 
critical coupling constant Acrit(-^) such that Hj\f is bounded below for all A < Acrit(-^) and 
such that 

mt — - — ^77^ — = — oo 

l/'eL2(R3iV) 
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for all A > Acrit(-^)- It was proven in [15] that, as N ^ oo, Acrit(-^) — ^ ^ait^ where A^j 
the critical constant for the Hartree energy 

Wee(v^) = / dx (1 - A)VV(.) dxd.M^)^^ . (8) 

J 2 7 \x-y\ 

In other words, A^^^ is such that <fHartree(9') > for all ip £ i?^/^(R^) if A < A^;^ while, if 

^ > KriV 

inf Suartrceif) = "OO . 

V5g//l/2(K3)_||^|| = l 

The criticality of the Hamiltonian ([7]) is a sign for the instability of boson stars when 
A = NG > A^j^. If the number of bosons in the star exceeds the critical value A'^crit = 
A^jj/G, the star collapses. 

Next, we focus on the properties of the evolution generated by the Hamiltonian 
defined in ([7]). For the subcritical regime, we show in [3], a joint work with A. Elgart, that 
the many body dynamics can be approximated (in the sense of Theorem I l.ip by the solution 
of the relativistic Hartree equation 

idt(pt = Vl-A(^t-A(T^* Iv^tM (ft ■ (9) 



Note that, in the subcritical regime, Lenzmann showed in |12] that Eq. ^ is globally 
well-posed in the energy space H^^'^{Mp). What about the supercritical regime? Since Hn 
is not bounded from below, it is not a priori clear how to define the one-parameter group of 
unitary transformations Uj\f{t) = e"*^'^* describing the time-evolution. To circumvent this 
problem, we introduce a tiny, A^-dependent, cutoff a{N) in the Coulomb potential, 

N A ^ 1 

ffS = E;^-sE^^7ri^^TlH^ (10) 

j = l Kj ' 

and we assume that a{N) — )• as A^ — )• oo. The regularized Hamiltonian H'^ is now bounded 
below for every A^. Therefore it can be extended (uniquely) to a self-adjoint operator on 
Hn, the unitary group UN{t) = e~*^^* is well defined and the Schrodinger equation 

idt'4'N,t = H^iljN^t (11) 

is globally well posed 0x1%^. On the other hand, since the cutoff vanishes in the limit 
A^ — )• 00, we may still expect the effective dynamics to be described by the semi-relativistic 
Hartree equation ([9]). 

It is important to notice that the criticality of the model can also be observed at the level 
of ([9]). For A > A^l^j^., the equation is still locally well-posed in H^f'^{M?) (for an arbitrary 
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initial data Lp in if^/^(M^) there exists a unique solution ipt in H^^'^ on the time interval 
t E {—T,T), for some T > 0). In general, however, the local solution cannot be extended 
to a global solution (i.e., one cannot take T = oc). In fact, for A > A^j^, it was proven 
by Frohlich and Lenzmann in [9] that there exist solutions of ^ which exhibit blow-up in 
finite time. This means that there are solutions ft of and < T < oo such that 



/2\ -^/^ 
dx (1 - Ay/^ft{x) \ oo (12) 



as t — )• T . Solutions of dH) exhibiting blow-up in finite time can be used to give a dynamical 
description of the phenomenon of gravitational collapse. 



The next two theorems from [16], a joint work with A. Michelangeli, show that, also in 
the supercritical regime, the solution ^pt to the relativistic Hartree equation ([9]) continues 
to approximate the many body dynamics until the time where ipt blows up (if ft does not 
exhibit blow up, then it stay close to the solution of the many-body Schrodinger equation 
on every finite time interval). 

The first theorem proves that, if ipt does not blow up in the time interval [—T, T], \(pt){ft\ 
is close to 7^j, as — )• oo, for all |t| < T. 

Theorem 2.1. Fix ip G H'^{M?) with \\ip\\ = 1 and T > such that 

K := sup ||(/Jt||^i/2 < oo (13) 
\t\<T 

where (ft is the solution of (0) with initial data (pt=o = Let ipN,t = e~*^^*93®^. Then 
there exists a constant C (depending only onT, \\ip\\H2, and k) such that 



Tr 



< (14) 



for all \t\ < T. 



The next theorem shows that if ft blows up at time T, then also the solution to the 
many body Schrodinger equation V'iv,* collapses as t approaches T, if, at the same time, 
N ^ oo sufficiently fast. This result justifies the use of the Hartree equation (jH) for the 
description of the gravitational collapse of boson stars. 

Theorem 2.2. Fix f G //^(M^) with \\ip\\ = 1. Suppose that > is the first time of 
blow-up for the solution ft of with initial data f (\\ft\\jji/2 < oo for all t G [0,Tc) and 
\\ft\\m/2 oo as t ^ T^). LetipN,t = e'^^^^V^^^ assume that, in [W\). a{N) > 
for some /3 > 0. Then, for every \t\ < Tc, there exists a constant Ct < oo such that 

11(1- A^ji/^V'iv.tll <C7t 

uniformly in N. Moreover, if N{t) G N for t G [0,Tc) is so that N{t) — t- oo sufficiently fast 
as t ^ Tc, we have 

||(1-A,JV4^^^^^^^||2^ r^(l_A)i/2^(;)^^^^^^ t^T-. (15) 
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To show these two theorems we use the techniques developed in |19j . These techniques 
were first introduced, in a shghtly different context, by Hepp in [10]. They are based on 
the representation of the many-body system on the Fock space, and on the use of coherent 
states as initial data. With respect to [19], the main novelty is that, to prove Theorem 12.21 
we need to show convergence not only in the trace norm (as in (I14p ). but also in the energy 
norm. In other words, we have to show that, as long as ipt does not blow up. 



Tr 



as — )• oo. In fact, as pointed out to us by R. Seiringer, the existence of < Tq < Tc with 
the property ||(1 — A^, J-^/^ VAr{t),tlP — )• oo as t — )• Tq" (with N{t) — )• oo as t — )• Tq~) follows 
from ()14p and from the semicontinuity of the kinetic energy (the kinetic energy of the limit 
is always smaller than the limit of the kinetic energy). However, it is only ()16p that allows 
us to conclude that Tq = Tc and therefore that the collapse of the many body wave function 
can really be described by the blow up of the solution of the Hartree equation ([9|). 



3 Dynamics of Bose-Einstein condensates 

Since the work of groups around Cornell and Wieman at the University of Colorado, and 
around Ketterle at MIT, see [U [2], Bose-Einstein condensation has become accessible to 
experiments. In these experiments dilute Bose gases are initially trapped by strong magnetic 
fields. Then, after cooling the gas to very low temperatures (of the order of nano-kelvin) , 
the traps are switched off and the evolution of the gas is observed. To understand these 
experiments it is important to find an accurate description of the macroscopic properties of 
the evolution of the condensate. 

The trapped Bose gas is described by the Hamiltonian 

N N 

hT = + V,^t{xj)) + N^V{N{x, - Xj)) (17) 

j=l i<j 



acting on the boson Hilbert space T-Ln = Tf(M ). Here T4xt is an external potential 
modeling the magnetic traps, V > is repulsive and of short range, and the interaction 
potential Vat (2;) = N'^V{Nx) scales with the number of particles so that its scattering 
length is of the order 1/A^. Recall that the scattering length of V is defined as 



Svrao 



j dxV{x)f{x) (18) 



where / is the solution of the zero-energy scattering equation 

A + ^f)/ = (19) 

with the boundary condition f[x) — )• 1 as |x| — t- 00 (it is then simple to check that /(x) ~ 
1 — (ao/|a;|) + 0(|x|~^) as |x| — ?■ 00). If / is a solution of (fT9]) . it follows by scaling that 
/Ar(x) = f{Nx) solves 

-A + iy^v) /7V = 0. (20) 
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This implies, in particular, that the scattering length of Vn equals uq/N. 

In [13], Lieb, Seiringer and Yngvason proved that, if -Etv denotes the ground state energy 
of Hn, En/N — )• min<^.||(p||=i i5gp(99) as — )• oo. Here 

£gp{^) = J dx (|Vv9(x)|2 + yext(x)|99(x)|2 + 47rao|(^(x)H 

for all if S L^(M^) is the so called Gross-Pitaevskii energy functional. In [13], Lieb and 
Seiringer showed that the ground state of H^"^^ exhibits complete condensation into the 
minimizer (pcp of £gPi in the sense that 

I'i^ ^\(I)GP){4>GP\ as N^oo, (21) 

where 7^'* is the one-particle density of the ground state of H^^'^^. 

What happens now when the traps are turned off? The system starts to evolve with 
respect to the translation invariant Hamitonian 

N N 

Hn = Y1 -^-j + N''y{N{x^ - xj)) . (22) 

It turns out that the macroscopic properties of the resulting evolution can also be ap- 
proximated by the same Gross-Pitaevskii which successfully describes the properties of the 
ground state. The following theorem is proven in [H O [7], a series of joint works with L. 
Erdos and H.-T. Yau. 

Theorem 3.1. Suppose that < V{x) < C(l + x^)^'^/^ for some a > 5. Let ipN G 
L^(]R^^) be a sequence of N -particle wave functions with \\4'n\\ = 1, so that, as N ^ 00, 
— ?■ \^){^\ for an arbitrary if G //"^(M^) (tp^ exhibits complete condensation) and so that 
[iI^m^Hn'^n) ^ CN (ipN hcLs finite energy per particle). Then, for every fixed t G M, the 
evolved wave function ipN,t = e~*^^*^7v still exhibits complete Bose-Einstein condensation, 
in the sense that 

iN]t^\^t){n\ (23) 
as N ^ 00. Here ipt is the solution of the Gross-Pitaevskii equation 

idt^pt = + S'naQ\ift?Lpt (24) 

with the initial data (pt=o = ^■ 

Observe that (j23p immediately implies that — )• \(pt) {iptl®^ for all > 1 (this is a 
general properties of densities converging to rank one projections). This theorem shows that 
the Gross-Pitaevskii equation can be used to describe the evolution of the condensates in 
the experiments we discussed above. In the rest of this section, I will explain the main ideas 
of the proof of Theorem 13.11 We follow the general strategy introduced in [20] to analyze 
the mean field limit; we need however to adapt the techniques to the present situation. 
The main difficulties are a consequence of the fact that (|22p does not describe a mean 
field regime, which is characterized by many weak collisions among the particles. On the 
contrary, the evolution generated by (I22p is characterized by very rare and very strong 
collisions (particles interact only when they are extremely close, at distances of order 1/N). 
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The starting point of our analysis is the study of the time evolution of the reduced 
density matrices associated with the solution of the Schrodinger equation tpN,t- It turns 

(k) 

out that the family Tnv ^ ~ l;---;-^) satisfies a hierarchy of coupled equations, 



commonly known as the BBGKY hierarchy: 



■a (fc) 



k k 

j = l i<j 
k 

+ iN- A;) ^Trfc+i N^ViN{xj - Xfc+i)),7S 



(25) 



t 



where Tr^+i denotes the partial trace over the degrees of freedom of the (A; + l)-th particle. 
From this hierarchy, we try to find equations for the limit points of the reduced densities, as 
— )• oo. Suppose, for example, that 7^^^ and 7^^^ are limit points of 7^^^ and, respectively, 

of 7^ J. Then, from ([25l) . with = 1, we may expect 7^^ and 7^\ to satisfy the equation 



idtlZ = [-^' r^lt] + Tr2 [boS{xi - X2),^Z 
where we used that, as — )• 00, 

{N -l)N^V{N{xi-X2)) ^bo6{xi-X2), with 60 = 



(26) 



V{x)dx. 



It turns out that (|26p is not correct. Taking the limit — t- 00 in (|25p . we ignored the 

(2) 

correlations developed by the two point reduced density 7]^ ^ on the length scale 1/A^. 
Using the solution /Ar(x) = f{Nx) of the zero energy scattering equation ([20]) to describe 

the correlations, we can try to approximate 7^^, 7^^^, for large but finite A^, by 



(1) 



7^^^(xi,X2;x'i,X2) ~ /Ar(xi -X2)fN{x[ - X2)-f^^'t{xi, X2; x[, X2) . 



(2) 



(27) 



Inserting this ansatz in (|25p (for k = 1), and using the definition (jlSp of the scattering length, 

(28) 



we easily find another (and, this time, correct) equation for the limit points 7,^\,7^\: 



loo.t 



+ 87raoTr2 



6{xi -X2),-fS,t 



Note that, /tv — ?• 1 as A^ — >• oo (in a weak sense). For this reason, /tv only plays an 
important role when it is coupled with the (very) singular potential V/v- Similarly to (|28p . 
starting from (|25p for k > 1 and considering the limit N ^ oo (taking into account the 
short scale correlation structure developed by 7^^^^), we conclude that an arbitrary limit 

point {7^\}fc>i of the sequence of reduced density matrices {'~f^j^\}^^i satisfies the infinite 
hierarchy of equations 



idtX 



OO.t 



n- 



-A 'Y^'^)' 



+ 87rao Trfc+i d{Q 



(k+i) 



■j - Xk+l), 700,1 



(29) 
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It is then useful to observe that the factorized densities 

tS = \n){^tf' (30) 
solve (f29]l if and only if ipt is a solution of the Gross-Pitaevskii equation (fM|l . 
The strategy to show Theorem 13.11 consists therefore of the following steps. 

1) Show the compactness of the sequence of families {7^{}^i with respect to an appro- 
priate weak topology. 

2) Prove that the reduced densities T^^^'', k>l, have a short scale correlation structure 
which can be described, in good approximation, by the solution /jv of the zero energy 
scattering equation ([20]) . 

3) Using the short scale structure developed by 7/^^^\ show that any limit point of the 
sequence {7^|}^i is a solution of the infinite hierarchy ([29|) . 

4) Prove the uniqueness of the solution of the infinite hierarchy ()29l) . 

Steps 1-4 conclude the proof of Theorem 13.11 because a compact sequence with at most one 
limit point must be convergent and therefore (j3U|) must be its limit. Details can be found 
in [H m m [7]. Recently, an alternative proof of Theorem 13.11 was proposed by Pickl in 
[TTJUS]. Pickl's approach also allows for the presence of (possibly time-dependent) external 
potentials in the Hamiltonian (j22p (but requires stronger conditions on the initial wave 
function) . 
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